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SUMMARY 


Transition of the Boundary layer fro® the laminar to 
the turbulent regime was investigated on the concave side 
of a plate with a radius of curvature of 2.!" feet* The 
critical Reynolds number was founc! to be considerably lower 
than on a flat plate and on the concave side cf a plate with 
a 2 O-foot radius of cu'-vat u^e previously investigated. It 
was furthermore found that, in agreement with the theoretical 

results of Gortier, R 9 J , here termed the “Gortlor param- 
eter, 11 is the proper critical parameter governing boundary 
leyer instability due to concave curvature. The critical 
parameter at transition was found to have a value of 9,0. 

The influence of pressure gradient and of an increased 
free— stream turbulence level on the position of the transition 
point on the concave side of the plate of radius of 

curvature has been studied. Small variations of the pressure 
gradient did not alter the value of the critical Gortler pa- 
rameter. This result is compared with similar measurements 
on the c cn v ex side of a plate of 2 0— f o o t radius of curvature. 
Increased tunnel turbulence lowered the value of the critical 
parameter at a turbulence level of 0.0 p percent to 6 C 0 

at a turbulence level of 0.3 percent. 

The investigation confirms the previous result that the 
mechanism of the breakdown of the laminar boundary layer is 
essentially different ori convex and concave boundaries, 

A discussion of the practical appl icabi 1 i ty of transi- 
tion measurements is gi ven 1 and the difference betweeen crit- 
ical Reynolds number corresponding to laminar instability and 
critical Reynolds number corresponding to transition is pointed 
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out r A definition of transition Reynolds number, based on 
the apparent shearing stress caused by the laminar oscilla- 
tions, i s given. In the case of flow along a flat plate, 
values ol transition Reynolds number are calculated approxi- 
mately for different intensities of the initial disturbance. 


INTRODUCTION 


..he larger pert of the California Institute of Technology 
investigation on the effect of curvature on boundary layer 
stability and transition has been previously reported ( refe r— 
ence l), -his previous research led to the conclusion that 
the mechanism of the breakdown of the laminar boundary layer 
is essentially different on concave and convex walls. The 
instability of the laminar boundary layer on convex boundaries 
was found to be brought about by plane disturbances, the so- 
called 1 o 1 lmi e n * Schlichting waves, w h i c h were first observed 
by Dryden, Schubauer. and $k rams tad on a flat plate. The 
i an g a of un s o able frequencies an d the am p 1 i f i c a t i o n c h a r act e r— 
is tics of the foilmie n— Sen lighting waves were found, within 
the experimental scatter, to be the same on the? converly 
curved boundaries as had earlier been observed by Schubauer 
and Skramstad (reference 2) on the fist plate. It was accord- 
ingly concluded that the transition point was essentially un- 
affected by convex curvature. These measurements extended 
over a range of effective curvature 0/r from 0 to ICT 3 ' 
where 0 is the momentum thickness of the boundary layer and 
r is the radius of curvature of the boundary, 

concave curvature was found to have a pronounced destabi- 
lizing eft ect on the laminar boundary layer, the critic a. 1 
Reynolds number decreasing with increasing effective curvature 
It was, therefore, concluded that the mechanism of transition 
on concave walls was different from that on flat an d convex 
boundaries. It was furthermore sr. own that these observations 
were in general agreement with theoretical results of Gortler 
(reference 3; which predicted a strong dynamic instability due 
to three-dimensional disturbances on concave boundaries. The 
measurements for the case of concave boundaries were however 
noc so extensive as for the convex case and covered only the 
range of effective curvature between 0 and 0,0001, One pur- 
pose of the investigations reported here is the extension of 
the previous measurements to values of effective curvature of 
about 0„ 001 on a concave boundary 


The influence of the pressure gradient along the 
boundary layer on the position of the transition point 


1 am in ar 
has long 
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been known. The measurements which have been carried out 
before were always confined to flat-plate flow or flow 
pact convex boundaries.. Based on Gortler's theoretical 
results It was suggested in the previous report that the 
influence of pressure gradient on transition should be less 
pronounced on concave boundaries. Consequently measure- 
ments of the effect of pressure gradient on the position of 
the transition point have been carried out and the results 
are presented in the present report. The effect of an in- 
creased free— stream turbulence level on transition was also 
studied. 

The investigations reported here confirm and supplement 
the results reported previously, and thus conclude the re- 
search program on the influence of curvature on boundary- 
layer transition carried out over a number of years at the 
California Institute of Technology under the sponsorship 
and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

The continuous interest, encouragement, and advice of 
Nr. Th. von Kansan and Dr. 0, B. Millikan is gratefully 
acknow ledged . The author also wishes to acknowledge many 
helpful discussions with Dr. C, C„ Lin, and the assistance 
of Mr, J, Laufer.. 


SYMBOLS 


x distance along surface of plate from leading edge 

y distance perpendicular from surface 

z distance along surface perpendicular from leading 

edge 

r radius of curvature of the boundary layer 


u 

mean velocity 

at a point 

in 

the boundary 

laye r 

U o 

mean velocity 

of the free 

stream 


u 

instant aneons 

x— c omponent 

of 

fluctuation 

velocity 

V 

instantaneous 

y— c component 

of 

f luc t uat i on 

velocity 

w 

ins t ant an eons 

z— c omponen t 

of 

f luctuat ion 

velocity 
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u ' 

V ! 
W ! 
U» 

U 


^ root— mean— square values of u, v f and w, respectively 

etc c , turbulence levex. usually expressed in percent 

d e n s i t y 
viscosity 


p = j± kinematic viscosity 
P 

q «= -ypUc^ dynamic pressure of the fr oe stream 


■t T 0 

h = y / — Bias ius ' nondimensi onal parameter 
v v x 

boundary layer thickness 


6 

6* 

e 

1 

r 

R 


boundary layer displacement thickness 
boundary layer momentum thickness 

effective curvature of the boundary layer 

Reynolds number 


R e > 


R 


x tr ; 


tr 


E e J~ 

t a = - p uv 
du 


Reynolds numbers besed on x ? 9 pnd 6*. 
respectively, P. x = U 0 x/ Vss x— Reynolds 
number, and so forth 

Reynolds numbers at the transition point 

OS 

Gor t le r parameter 
apparent shear 


t t . - ti- 


dy 


laminar shear 
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k = correlation coefficient 

u * v > 

Cf L laminar sk in— f r i c t i on coefficient 

APPARATUS AND METHODS 
Wind Tunnel and Test Section 


The general layout of the wind tunnel is shown in 
figure 1. A detailed description can be found in reference 
1. The measurements were carried out on the concave side of 
a smooth glass nlete of 2.f-foot radius of curvature. This 
plate was set in the center of the high (.* ft) and narrow 
(7.F in.) test section shown in figure 2. The section con- 
sists essentially of two plate— glass plates forming cylinders 
concentric with the center sheet , on which the measurements 
are made. These plates are mounted in wooden frames which, 
by means of a screw adjustment, allow a variation of the 
width of the channel and of the angle of atte.ck of the center 
plate„ This adjustment makes it possible to alter the pres- 
sure gradient along the test section. Owing to the geometry 
of the section it is, however, difficult to set for a pre- 
determined simule pressure distribution, for example, a 
linear distribution. The investigations on the effect of 
pressure gradient on transition were, therefore, rather tedious 
and the investigated range, of pressure gradients limited. 


Free— Stream Turbulence Level 

The normal free— stream turbulence level of the tunnel 
was the same as in the previous measurements. The values 
wer e 

u i _ . 

— = 0, OF percent 


— — = — — = 0.12 percent 

U 0 Uo 

The turbulence level in the test section is mainly controlled 
by the precision screen "S" (fig. 1) in the pressure chamber. 
To investigate the position of the transition point at higher 
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turbulence levels of the free stream, the tunnel turbulence 
was artificially raised in the following way: Strips of 

celluloid tape about 1 inch wide and F inches long were 
glued at regular intervals of about 3 inches on the down- 
stream side of the screen S. These strips produced a 
uniform turbulence level of 

u • _ _ 

— =• 0,3 percent 

u 0 

in the test section. By removing every second strip, and so 
forth, the turbulence level could then be lowered to the 
normal free stream level. 


Traversing Me chon ism 

The traversing mechanism used in this wind tunnel is 
described in detail in reference 1. The motion of the 
instrument , for example, hot-wire anemometer, is continuous 
in the x and y directions and is remotely controlled. 


H o t— ST ire An e m 022 eter 


The mean— speed distribution in the boundary layer was 
measured with platinum wires of 0, 0005— inch diameter and 
about 5—mi llime t er length. The turbulence level and the 
position of the transition point were determined with 
0 n 00024— inch thick and about 2— millimeter long wires. 

Mean speed was always measured using the constant-resistance 
method „ For measurements of velocity fluctuations the wire 
was calibrated using the electrical oscillator method (refer 
ence 4 ) and the amplifier properly compensated by means of 
an inductance-resistance circuit. The amplifier response 
is flat between less than F cycles up to 8000 cycles, 


Measurement of the Pressure Distribution 

The pressure distribution along the test section was 
measured by means of a small static tube mounted on the 
instrument carriage at a distance of about 1 centimeter 
from the test plate. Thus, the pressure gradient close to 
the edge of the ooundary layer was measured. The irregular 
waves' in the distribution (fig o 3) rare due to local change 
in the radius of curvature of the plate. 
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Determination of the Transition Point 


The transition point was determined by means of the 
hot-wire anemometer. The hot wire was mounted on the 
instrument carriage a few thousandths of an inch from the 
surface of the plate. The carriage was then moved slowly 
in the direction of the mean flow from the leading edge 
downstream; the first appearance of the turbulent bursts 
(see reference l) on the oscilloscope screen was taken as 
indication of transition. 


The test— section walls were adjusted to give a pres- 
sure gradient as close to zero as possible* The best re- 
sult which could be obtained is shown as "A" in figure 3. 

The pressure gradient is seen to be close to zero over the 
first 90 centimeters of the test plate. Since transition 
was found to occur at distances from the leading edge of 
less than £0 centimeters* this pressure gradient A (fig. 3) 
could, w/ith sufficient approximation, be taken as 11 zero" 
gradient, Tnr velocity distribution in the laminar boundary 
layer was then measured. The result is shown in figure 4. 

It is seen, that the profile follows s rather closely, the 
Bias i us flat--plate profile. There exists a slight system- 
atic deviation from the Blasius profile especially in the 

range between p = 3 end p ~ 6 , The measured values of — — 


are slightly lower than the ones given by the Blasius dis- 
tribution. 

The main purpose of the measurements was to establish 
whether computations of the momentum thickness 9, the 
Reynolds number based on 9, and so forth could be computed 
from the Blasius solution. For computation of these data 
applied to transition measurements, the slight deviation 
seen in figure 4 is imma,ter ial and consequently the values 
of 9, 8* , and so forth, used later on, have been computed 

from the Blasius formula, that is, 


MEASUREMENTS AND RESULTS 


Mean Speed Dis tribat ion 
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Effect of Curvature on the Position 
of the Transition Point 


transition was measured with the hot — wire anemometer. 

By varying the velocity of the free stream, the range of 

effective curvature from -0.F x l(T 4 to -1.1 x 10~ 4 could 
he investigated. The results of the tests are shown in fig- 
ure F where the critical Peynolds number H x tr is plotted 
versus It is seen that R^. fells from about 4 x 10 5 

at i = _ 5 x 1C -4 to l.F x 10 5 at - = - 11 x 10~ 4 
r r 

This result agrees very well with previous measurements 
( reference 1). Figure 6 shows the effect of curvature on 
transition for effective curvature ranging from -0. 001 (con- 
cave)) to +0.001 (convex). Here R 0 tr is plotted versus 
0 

Jt is seen that the decrease of the critical number found 

on the concave side of the 20— foot plate continues with higher 
effective curvature. 


According to Gortler (reference 3) R 0 /- is the char- 

V * 

acteristic parameter for the three-dimensional instability 
on concave walls. I? the value of this parameter at transi- 
tion is plotted versus effective curvature (fig. 7a), it is 
seen that within the rather large scatter the value of 


V 




0 

does not systematically vary with — , The average 


f / 0\ 

value of /R 0 / - J found is about 9,0. The value of the 
^ / 1 r 

parameter as measured here is somewhat larger than the value 
of v.3 found previously on the concave side of the plate 
of 20— foot radius. There exists a possible explanation for 
this difference in the value of the Gortler parameter at tran- 
sition which can best be seen from a plot such as figure 9. 


H 

Instead of plotting the Gortler parameter as a function 
of effective curvature P x tr is plotted as n function of 

-g. Since R x ~R 0 ' C ) R 0 j~~&- - constant corresponds to 

v r 

straight lines through the origin in this plot. The measured 
values of t r ^ or the two test sections of 20-foot and 

2.5— foot, radius are presented in figure 8, and in addition, 
the line R 0 = 9.4 x 10^ corresponding to transition on a 
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flat plate in the wind tunnel is shown. Assume that transit 

M 

t ion due to Gortler vortices occurs for all concave curva— 

/ a 

tures at the value of J — found on the plate of 2 . 5 * — 


foot radius. All transition measurements should thus fall 


on the straight line 


.... /r> 

V y J ry 


R 


= 9.0 (fig. 8). Trans i- 


t r 


tion on the flat plate, that is, transition caused by the 
two-dimensional perturbation, was found to occur at 

S 0tr = X Since the Tollmien— Schlicht ing wave is 


but little affected by smell curvature, it can be expected 
that transition caused by this di sturbance ~on slightly con- 
cave plates occurs at about the same Eg . Hence, starting 
from the flat plate, all transition measurements should 

follow the straight line R 0 = 9.4 x IQ 2 in figure 8. At 
~ = - 11 x 10 3 the line ( /J- ) =9.0 and R Q = 9.4 x 

0 \ v 7 r / t r 9 

10""* intersect; that is, at effective curvatures larger 

than —11 x 1 0 3 transition is caused by the three-dimensional 
disturbance. At smaller effective curvature the Tollmien— 
Schlichting wave is more unstable and thus is responsible for 
transition. It appears that this is not the case, but that 
a jnore or less continuous change from the transition due to 
Gortler vortices to transition caused by Tollmien— Schl icht ing 
waves takes place. Since the laminar boundary layer on con- 
cave walls is unstable with respect to small perturbations 
of both kinds, such a continuous change is possible and 
probable. types of perturbation are excited by the 

external free— stream turbulence, traveling downstream, that 
is, toward higher values of Eg both the waves or the vor- 
tices of a given wavelength increase in amplitude when they 
pass through their respective instability region in the 
Rg*_ a diagram (see reference l). A transfer of energy from 

the one type of perturbation to the other is therefore quite 
possible and such a transfer becomes unimportant only if one 
of these perturbations is considerably more unstable than the 
other, that is, in the limiting cases of small and large ef- 
fective curvature, A more exact analysis of such an energy 
transfer is impossible without a detailed knowledge of the 
amplification of both types of perturbation and furthermore 
of the excitation due to free stream turbulence. 
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Effect of the Pressure Gradient on Transition 

It is well known that a pressure gradient along a 
laminar boundary layer strongly influences the position 
of the transition point* The so-called lo*w— drag or laminar 
flow airfoil represents an application of this result, All 
investigations on the influence of the pressure gradient 
have previously, however, been confined to measurements on 
flat or convexly curved boundaries, that is, to boundary 
layer flows, where the Tollmien— Schlichting in s tabi 1 i ty is 
predominant. Based on the theoretical analysis of G-or tier 
(reference 3) it was suggested previously (reference l) that 
the influence of the pressure gradient on transition should 
be less pronounced on concave walls where the boundary— layer 
instability should be due to the three-dimensional vortices. 

To investigate this question, measurements of the in- 
fluence of pressure gradient cn transition were carried out 
on the convex side of the plate of 20— foot radius and on the 
concave side of the plate of 2 foot radius. Thus, a typi- 
cal Tollmien— Schl ichting and Gortler case could be compared. 

The pressure distributions for which transition was investi- 
gated are shown in figures 3 and S 9 respectively. Before a 
discussion and comparison of the measurements can be presented, 
the question of the proper parameter for the pressure gradient 
and the proper critical number for transition must be dis- 
cussed* 


It is sometimes useful to present transition data in 
the form of a critical Reynolds number R x based on the 
distance from the leading edge, that is, from the stagnation 
point. However, R x certainly is not the proper physical 
parameter, and it is preferable to base the critical number 
on some measure of the local boundary layer thickness* The 
momentum thickness appears as the best choice since it is 
closely connected with the shearing stress at the wall and 
thus with the slope of the velocity profile at the w a 1 1 f In 
this connection the instability of a velocity "profile with 
respect to To 1 Xdj i en — Sc hi i c ht i ng waves is known to depend very 
much on the slope (and curvature) of the profile in the neigh- 
borhood of the wall, Gortler found that for the thrse-di;iea- 
sicnal disturbance also momentum thickness is the most suitable 


length 
to plot 

and for 


for the critical parameter. Thus 
for flat and convex surfaces $he 

the concave, surfaces ( P. A /l ) 

V 0 V r /tr 


effective pressure 

that is, — AP. i n 
q ax 


gradient* To represen 
dimensionless form, a 


it appears reasonable 
critical number E $ r 

as functions of some 

t the pressure gradient, 
length is again needed. 
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The effective force acting on an element of the boundary 

} 

layer is the force due to the pressure gradient. ~ . 

q dx 

The boundary layer thickness 8 entering here has core 
kinematic tha.n dynamic character and thus the displacement 
thickness 8* appears as the logical choice. This can 
also be seen from \ on. Kerman's momentum equation of the 
boundary layer 


T o _ 1 d ( 9 a ) _ dp 
q q dx q dx 


where also the conciliation 


5 * dp 

— — ^ appears, 

q dx 


The results of the teasureaents are, therefore, renre— 

sented in the form R^tr ^ » respectively, 

6* dt> V v r ytr 

versus — — — . It is. of course, not to be expected that a 

q d x 

complicated phenomenon, such as the laminar instability for 
flow with piessure gradient, can be described by one critical 
parameter alone, Under the influence of the pressure gradient, 
the velocity profile of the layer will, in fact, change contin- 
uously from the stagnation point up to eventual separation. 

The instability zones for snail perturbations and the ampli- 
fication characteristics thus vary continuously downstream* 
Transition occurs if pert urbat* i ons are amplified to a suffi- 
cient degree. Since the total amplification is thus an inte- 
gral effect depending on all instability zones passed by the 
perturbation on the way from its origin downstream to the 
transition point, it seems hardly possible to give a repre- 
sentation involving only one critical number such as R^ m 
Only where the instability zones are little affected by a 
change in the profile an, for examule, according to Gortler 
(reference 3), in the case of flow past concave boundaries 
can it be expected that it v/ill be possible to represent the 
instability character by only one critical parameter, in this 


case Hq / 

° v 


r ' 


Except for pressure gradient ?, C n on the convex 


side of the plate of 2 0— fort radius, the momentum thickness and 
d isplac ement thickness were again computed using the Blasius sc 
lution. The influence of these comparatively s~ail gradients o 
the bound arv—i^yer thickness was found to be negligible within 
the scatter of the transition measurement s . In the case of gra 
* lent C in figure 9 f the momentum thickness end displacement- 
thickness were computed from the solution of the laminar 
boundary— layer equations for the case of a linear decrease 
of the free— stream velocity. The experimentally determined 
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pressure distribution was approximated by a parabola, and 
the momentum and displacement thicknesses were determined 
from the tabulated solution (reference F ) ,, It is seen from 
figure 10 that for the case of convex curvature E 0 tr de- 


creases with positive and increases with negative pressure 
gradient as expected. The largest decrease with positive 
pressure gradients occurs close to zero pressure gradient. 
Whether or not this is true for negative gradients remains 
undecided since no measurements at larger negative values ©f 
iLUL ikE could be carried out due to the limited Reynolds number 


q dx 

of the test plates. The rapid change of Rq ^ r in the neigh- 
borhood of zero pressure gradient and the comparatively small 
change at larger positive values of — is interesting. 

Thus it is indicated that the Blasius solution is singular 
in the sense that a slight variation of the profile from the 
Blasius solution does cause large changes in the stability 
character. This is not unreasonable since the Blasius solu— 

" ti 

tion is indeed a singular solution 'with — — 5 = 0 at y = 0 . 

dy* 

A theoretical analysis of the instability character of pro- 
files where the pressure gradient is not zero but is small 
will be required for a physical explanation of this pressure- 
gradient effect. 


The absolute values given in f igur e 10 for P-Q^ r as a 

function of — — -•£ are certainly not universally applicable 
q dx 

since they are strongly influenced by the wind-tunnel charac— 
acteristics, that is, turbulence level, turbulence spectrum, 
and so forth (see also reference l). The relative variation 
of the critical Reynolds number at transition with pressure 
gradient, however, is believed to be quite general. Conse- 
quently, it is to be expected that in practical applications, 
such as the laminar flow wing, and so forth, the movement of 
the transition point when the pressure gradient is altered 
is most pronounced in the neighborhood of zero gradient. 


Figure 11 shows the effect of pressure gradient on 
transition on the concave side of the plate of 2 . F— foot 
radius. It is seen that the critical parameter for transi- 
tion of the boundary layer decreases with increasing pres- 
sure gradient. This effect hardly exceeds the scatter of 
the points, however, and is not believed, to be actually 
systematic. A comparison of the values of the critical pa- 
rameter for transition at positive gradients with the range 
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of values at zero pressure gradient shows that the apparent 
decrease at positive pressure gradient lies within the ex- 
perimental error. The comparatively large scatter of the 
points is due to the inherent difficulty of transition meas- 
urements on a plate with a small radius of curvature and 
also to errors in the determination of the boundary layer 

thickness which enters both parameters Eg f ~ and — 

The measurements were taken at small values of R x because 
of early transition, and thus 9 and 5* vary comparatively 
rapidly with x. An error in x tr thus produces larger 
errors in 0tr and ®*tr then appears in the case of the 

convex plate where ^x^ r was five to ten times larger. The 
present investigation accordingly indicates that the influenc 
of pressure gradient on transition at concave surfaces is 
either zero or at most very small. 


Effect of Free Stream Turbulence on Transition 

That free— stream turbulence has a strong influence on 
the position of the transition point is well established. 

The first systematic investigation of this effect was made 
by Schubauer and Skramstad (reference 2), who measured the 
effect of varied free stream turbulence level on transition 
of the boundary layer on a flat plate. 

In the present investigation, transition was measured 
on the concave side of the plate of 2. F— foot radius for 
three turbulence :.le vels of the free stream. The results 
are given in figure 7. It is seen that an increased turbu- 
lence level leads to earlier transition. The critical number 

(jSg J decreases from 9.0 with ~~ = 0.06 percent to 

6.0 with = 0,3 percent. The magnitude of the turbulence 

u o 

effect is of the same order as in the case of a flat plate 
(reference 2). Thus it appears that the three-dimensional 
disturbances, like the Tollmien— Schli cht ing waves, are ex- 
cited by turbulence in the free stream. A more detailed 
discussion of the turbulence effect meets, as in the case 
of the Tollmien— Schl i cht ing waves, with the difficulty that 
the phenomenon depends on the spectrum (or on the scale) of 
tunnel turbulence, and furthermore on the amplification char- 
acteristics of the perturbation along the path traveled by 
the perturbation. 
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LAMINAE INSTABILITY AND TRANSITION 


Tha preservation of laminar flow up to very large 
•Reynolds numbers is of considerable practical interest. 
Experimental investigations of transition from the laminar 
to the turbulent regime have, therefore, been undertaken 
with a double purpose: To give an understanding of the 

mechanism of transition or at least of the breakdown of 
1 am inar motion, and. to permit prediction c f the occurrence 
of transition for a given practical case such as a certain 
airfoil. The first part of this problem has found a com- 
plete solution: Recent experiments (leference 1 and espe-r* 
dally reference 2) confirmed, in general, the results of 
earlier theoretical 1 nv as t i gat i ons of Heisenberg, Tcllmien, 
Schlicht ing ? G-ort ler ? an d others. On the basis of this e x— 
peri mental and theoretical evidence, and also of earlier 
work of Taylor on the flow between rotating cylinders, there 
appears to be no doubt that laminar flow always becomes un- 
stable for a sufficiently large Reynolds number. This criti- 
cal Reynolds number will depend on the form of the laminar 
motion, for example, laminar boundary -lay er flow, laminar 
pipe flow, and so forth, and on the type of disturbance, for 
example, two-dimensional waves or three-dimensional vortices, 
Inst, anility in the sense used here means that for all Hej n— 
olds numbers above the critical one there will exist a ran^e 
of, say : wavelengths or frequencies such that a disturbance 
of the proper frequency or wavelength will increase in 
ampl i tude . 

Thus the mechanism of laminar instability is clear and 
the critical Reynolds number thus defined can be predicted* 
However, for practical purposes, it is desirable to predict 
a different Reynolds number’ namely, the Reynolds number at 
wh ich the flow becomes turbulent. There is a. distinct di f — 
ference between these two Reynolds numbers; if the first 
Reynolds number R 1 , that is, the critical number in the 
sense of the small perturbation theory, is reached, amplifi- 
cation of certain disturbances begins. The second, or ,l prac— 
tical 11 critical Reynolds number R^ means that amplification 
of disturbances has already taken place and to such an extent 
that complete breakdown cf the laminar motion occurs. 

The critical Reynolds number R 1 is, in general, so low 
and the difference between R 1 and so large that for 

practical applications the prediction of R 1 only is of 
little use. For example, the laminar boundary layer along 
a flat plate in the absence of a pressure gradient becomes 
theoretically unstable at a Reynolds number based on the 
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distance from the leading edge, R = 6 x 10 4 . Experiments 
in low- turbulence tunnels (references 1 and 2) confirmed this 
number but showed also that actual transition occurred only 
at Reynolds numbers from 2 X 10 6 to 2.8 X 10 6 . Hence, the 
larger part of the boundary layer is in most practical cases 
in an unstable state, that is, in a condition similar to that 
of a supercooled liquid or a supersaturated vapor. The Rey- 
nolds number R 1 corresponds in this analogy to the melting 
^r condensation temperature T 1 . The Reynolds number R 2 
corresponds to the temperature T 2 at which under given ex- 
perimental conditions the substance a.ctually solidifies or 
condenses. Consequently R i is a definite number for a cer- 
tain type of laminar motion, cor respondi ng to the definite 
melting temperature of a certain substance. On the other 
hand, R s depends not only upon the type of laminar motion 
but also upon the intial disturbances present in the laminar 
layer and therefore on the experimental setup, for example, 
the free stream turbulence level. The a mp lification of these 
initial disturbances in the unstable region can be computed 
from the small perturbation analysis, provided the magnitude 
of these disturbances does not exceed the range of this line- 
arized theory* Let a(R x ) denote the ratio of the amplitude 
of a given oscillation at a Reynolds number R x to that of 
the same oscillation .iust upstream of the instability zone. 

The value a(R x ) depends, of course, on the frequency of the 
oscillation (see, for example, references 1 and 2). It is now 
possible to compare a(R ) for various frequencies and obtain 

the maximum possible amplification a(R x ) at a given Rey- 

m a x 

nolds number R x . Such a c omnut at i o n whs carried out by 
Schlicht ing (reference 1) for flat-plate flow. Schlichting 
found in this case 

= °.55 .‘—x. 0 -*, (D 

Comparing his result with measurements of Hansen and Gebers, 
Schlichting noted that at transition the ,( most dangerous” 
frequencies had been amplified four to nine times. However, 
if the values of R x found in recent experimental work are 
inserted into equation (1), a(R x ) becomes of the order of 
10 5 . This point was emphasized by a. I. Taylor in IQ38 
(reference 8) a.nd at that time was considered as strong evi- 
dence the validity of the small-perturbation theory. 

This apparent difficulty in applying equation (l) is mainly 
caused by the lack of a clear definition of R p , the Reynolds 

number of transition; thus, in ^rder to bridge the gap between 
R i it is first necessary to define R 2 in such a way 

that it can be related to the small perturbation theory and 
hence to R 1 t 
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The following definition appears to be reasonable: 

"S j is the Reynolds number at which the apparent shear 

T A = — p uv due to amplified boundary layer oscillations 
at my point in the boundary layer becomes equal to the 


laminar shear 
with H, 


T L = p 4“ in the boundary layer," Beginning 
'“*■ • y 

^ , certain disturbances will be amplified. The per- 
turbation components u and v are correlated and con- 
sequently give rise to an apparent shear. The magnitude of 
this apparent shear increases as the amplitude of the oscil- 
lations increases. When this apparent shear becomes of the 
same order as the laminar shear of the undisturbed mean mo- 
tion the mean profile can evidently no longer remain unaltered. 
Hence. R^ defined in this way can be expected to be reason- 
ably close to the Reynolds number at which the osc i 1 lat i on has 
an effect on the mean flow sufficient to modify radically the 
mean velocity profile. 


The shear which exists in the undisturbed laminar motion 
can be computed as a function of R x for any given case. 

The apparent shear caused by amplified boundary layer oscil- 
lations T^, is more difficult to compute since depends 

upon the initial disturbance in the laminar layer and the 
ampl i fication characteristics. 


It is , however, possible to bring into a fs>rm 

where its general behavior as a function of R x becomes 

more evident and where in certain cases, at least the order 
of magnitude of and thus of can be obtained by 

crude approximations. 


If the correlation coefficient k = 

! i u T v 1 

fluctuation levels and are introduced 

U,> 


an d the 

becomes 


‘A 


uv = 


- P 




111 


v 1 
U n 


u. 


It is furthermore convenient to introduce a factor b by 
writing 


u t _ v_» 


and thus 


NACA ACR IT o . 4,728 


17 




( 2 ) 


Using the amplification factor 

. . . /niN 

to its initial value, ( — I 


a(R), £L can be related 




just upstream of the ampli— 


f icati on sons: 



f u i 

2 k b ( — 



s 


(3) 


Equation (3) is still quite general, will, of course, 

vary across the boundary layer and its magnitude will also 
depend upon the oscillation frequencies and on the Reynolds 
number, Hovrever, only the max im um value which can 

reach at a certain Reynolds number is necessary for the 
present discussion, since an attempt was made to find the 
Reynolds number at which r A first becomes equal to T L . 

C onse quent ly , 


as function of the Reynolds number R Xo Now, the distri- 
bution of t A across the boundary layer has a maximum close 

to the so-called "critical layer" which is usually near the 
wall. Eence, for the transition criterion T L shall be 

replaced, approximately, by its value at the wall. Thus, 
in terms of the usual laminar skin friction coefficient, 
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Ifl 


A max h _ 2 

t L c f L 


(k b * [aU)A 

l \Pj 1 J J 


(P) 


max 


R ... is then determined by 


•i k b 

i 


(f) 

\ U o/i 


— 1 


( 6 ) 


max 


Application to Flat Plate Plow 

The results obtained by Schlichting (references 6 and 7) 
can now be used to evaluate equation (F) or (b) numerically, 

Equation (l)* gives a(R) . The correlation coefficient 

max 

k has been computed by Schlichting (reference 7) at one 
point on the lower branch and one point on the upper branch 
of the instability limit curve. Comparing the results at 

these two points it appears that k. differs but little 

max 

and thus the variation of k with Reynolds number is 

apparently very small. Since Cf and especially a(R x ) 

L 

vary rapidly with R x 1 it is sufficient to replace k(R) 

in equations (f) and (b) by a constant value. The maximum 
value k == — 0.18 found by Schlichting is thus used in equa- 
tion ( F ) , The same argument apparently holds for b$ in 
the neighborhood of the critical layer, b is of the order 
0.1. This value also agrees with measurements (reference l). 
The skin friction coefficient for the Plasius layer is given 

by 

Cf L = 0.664 

Thus equation (?) finally becomes: 


*The use of Schi icht ing ' s amplification function in the 
present analysis is not strictly correct; the oscillations are, 
of course, no longer "small perturbations" in the sense of the 
linear theory if becomes of the order of T The error 

involved is, however, probably not very large and the order cf 
magnitude of fU, which is of main concern here, will be un- 
affected. 


NACA ACR No, 4J28 


19 


V - H - o. oie (" 




& 7L e 


i • n ^ i o 


(7) 


Figure 12 eives 


A m a x 

tl 


plotted against R x for different 


values of the initial fluctuation level. In considering 
the magnitude of the initial disturbance, it should be kept 

in mind that ( - u 1 N - 


Vu 


is the fluctuation level with the most 


o 


is only- 


fraction of the total fluctuation level in 


dangerous, omit is, most highly amplified, frequency. Thus 

(si) 

V,/i 

the laminar layer upstream of the amplification zone. 

Two facts are evident from figure 12. 

T 


(l) The slopes of the curves versus 



Hence R, does not depend too much upon the exact 
t 

value of —^—2—— which is chosen to define transition. 
T L 


R x even with considerable ini— 


( 2 ) Quant i ty R £ > 
tial fluctuation levels. (For example, w i t h 

1 percent, R z = 10 R 1 . ) 


(s!) 
w i 


The next step would logically involve relating the 
quantities entering equation (f) to factors known to influ- 
ence transition. Especially the relation between 

and the free stream turbulence level is very im— 
i 

portant , 

At the present time only a few general c on s id er at i ons 
can be given in this connection and the main effects classi- 
fied, It is hoped to obtain more complete relations in the 
future. 

( 1) Eff e ct of free— stream turbulence Turbulence of the 



NACA ACE No. 4J28 


20 


free stream will affect the quantity 


o 


effect will depend, upon the free stream level 


only. The 

Vu 0 > f 


and upon the frequency distribution or upon the scale of 
turbulence. The relation between ) fiod 

vu 0 ; f 

not simple and requires further study. 

(2 ) Effect, of pressur e g ra di ent.— The pressure gradient 
will affect at least two quantities in equation (3)j namely, 
Cf T and the amplification function a(H). An adverse, that 

is, positive gradient causes a more rapid decrease of Cf^ 
with Reynolds number to Cf T = 0 at separation, and as 

Schubauer and Skramstadt (reference 2) found, a larger 
amplification a(R), 


(3^ Effect of _r oughnes s . - Small roughness elements, 

/u ' \ 

like free stream turbulence, will mainly affect <— / . 

vu o'i 

The value of (— ) corresponding to a given roughness 

^ U o y i 


element will depend upon the element's height and also upon 
the width, since the latter influences the frequency distri- 
bution of the disturbance caused by the presence of the 


element. 


C ONC LUD I NG RE KAR K S 


The results of measurements of boundary layer transi- 
tion on the concave side of a plate with a 2. F— foot radius 
of curvature confirm the previous result that the mechanism 
leading to the breakdown of laminar flow differs essentially 
at convex and concave boundaries. Transition on convexiy 
curved surfaces is due to the two-dimensional Tollmien— 
Schlichting type of instability as is the case with a flat 
plate. Transition on concave boundaries appears to result 
from the three-dimensional Gortler type of instability, and 
is thus related to the instability occurring in the flow 
between rotating cylinders. 
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No influence of pressure gradient on transition on the 
concave wall was found. The influence of a positive prsssure- 
gradient on transition on convex boundaries is most pronounced 
in the neighborhood of zero pressure gradient,. Tree— st ream 
turbulence appears to have about the same effect on the posi- 
tion of the transition point in concave and convex boundary 
layer flew. 


Tne mean velocity profiles in the boundary layer of a 
concave plate with an effective radius of curvature around 

ICT 3 was found to be very close to tne Blasius flat— plate 
d i s t r i but i on . 

A discussion of the breakdown mechanism of laminar flow 
shows the difference between the Reynolds number E., corre- 
sponding to the beginning of .instability, in the cense of a 
small perturbation theory, and the Reynolds number R. 
corresponding to transition. The former number is well 
defined and can be predicted for a given case with good 
accuracy, The latter can be defined by the condition that 
at R the maximum apparent shear due to the laminar bound- 
ary layer oscillations becomes equal to the laminar shear- 
ing stress at the wall. for the undisturbed mean motion, A 
rough estimate of R £ fer the case of flat plate flow shows 
that even with considerable initial disturbances is 

much larger than R 1# 

California Institute of Technology, 

Pasadena. Calif., October 28, 1944, 
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Figure 1.- Flow diagram of the boundary layer research tunnel. 
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Figure 3.- Pressure gradients along the concave side of the plate with 2.5-foot radius 
of curvature for which transition was investigated. 
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Figure 5.- Reynolds numoer of transition against the 
side of tho plate witn Z. 5-foot radius of 


effective curvature on the concave 
curvature . 
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Figure 6.- Effect of curvature on transition. The arrows point to the value of Rq 

found for the flat plate. tT 
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Figure 7.- Transition measurements on the concave aide of the plate with 2.5-foot radiue of 
curvature at three different turbulence levels of the free stream. 



Figure 8.- Transition measurements on the concave sides of the plates of 2 5 and 20-foot 
radius of curvature. 
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Figs. 9,10 



Figure 9.- Pressure gradients on tne convex aide of tne plate of 20-foot radius for wnicn 
transition was investigated. 
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Figure 10.- Effect of tne pressure gradient on traneition for flow past the convex side 
of tne plate with 20-foot radius of curvature. 
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Figure 11.- Effect of the pressure gradient on transition on the concave side of the plate 
witn 2.5-foot radius of curvature. 





